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Abstract. We address the question of a Berry-Esseen type theo- 
rem for the speed of convergence in a multivariate free central limit 
theorem. For this, we estimate the difference between the operator- 
valued Cauchy transforms of the normalized partial sums in an 
operator-valued free central limit theorem and the Cauchy trans- 
form of the limiting operator-valued semicircular element. Since 
we have to deal with in general non-self-adjoint operators, we in- 
troduce the notion of matrix-valued resolvent sets and study the 
behavior of Cauchy transforms on them. 

1. Introduction 

In classical probability theory the famous Berry-Esseen theorem gives 
a quantitative statement about the order of convergence in the central 
limit theorem. It states in its simplest version: If (Xj)jgN is a sequence 
of independent and identically distributed random variables with mean 
and variance 1, then the distance between S'„ := -^^{Xi -|- ■ ■ ■ -|- X^) 
and a normal variable 7 of mean and variance 1 can be estimated in 
terms of the Kolmogorov distance A by 

'n 



where C is a constant and p is the absolute third moment of the vari- 
ables Xi. The question for a free analogue of the Berry-Esseen estimate 
in the case of one random variable was answered by Christyakov and 
Gotze in j2] (and independently, under the more restrictive assumption 
of compact support of the Xj, by Kargin [lOj): If {Xi)i^^ is a sequence 
of free and identically distributed variables with mean and variance 1, 
then the distance between Sn '■= -^{Xi -|- ■ ■ ■ -|- Xn) and a semicircular 
variable s of mean and variance 1 can be estimated as 

A/c \^ l"^3| + ^AnI 



n ' 
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where c > is an absolute constant and and 7714 are the third and 
fourth moment, respectively, of the Xi. 

In this paper we want to present an approach to a multivariate ver- 
sion of a free Berry- Esseen theorem. The general idea is the following: 
Since there is up to now no suitable replacement of the Kolmorgorov 
metric in the multivariate case, we will, in order to describe the speed 
of convergence of a (i-tuple {Sn \ . . . , S^^) of partial sums to the limit- 
ing semicircular family (si, . . . , s^), consider the speed of convergence 

oip{Sn \ . . . , Sn^) to p{si, . . . ,Sd) for any self-adjoint polynomial pin d 
non-commuting variables. By using the linearization trick of Haagerup 
and Thorbj0rnsen |51 E], we can reformulate this in an operator- valued 
setting, where we will state an operator-valued free Berry-Esseen theo- 
rem. Because estimates for the difference between scalar- valued Cauchy 
transforms translate by results of Bai JT\ to estimates with respect to 
the Kolmogorov distance, it is convenient to describe the speed of con- 
vergence in terms of Cauchy transforms. On the level of deriving equa- 
tions for the (operator-valued) Cauchy transforms we can follow ideas 
which are used for dealing with speed of convergence questions for ran- 
dom matrices; here we are inspired in particular by the work of Gotze 
and Tikhomirov |4j, but see also [IJ. 

Since the transition from the multivariate to the operator- valued set- 
ting leads to operators which are, even if we start from selfadjoint poly- 
nomials p, in general not self-adjoint, we have to deal with (operator- 
valued) Cauchy transforms defined on domains different from the usual 
ones. Since most of the analytic tools fail in this generality, we have to 
develop them along the way. 

As a first step in this direction, the present paper (which is based on 
the unpublished preprint [14j) leads finally to the proof of the following 
theorem: 

Theorem 1.1. Let (C,t) be a non- commutative C* -probability space 
with T faithful and put A := Mm(C) ® C and E := id®T. Let (X)jgN 
be a sequence of non-zero elements in the operator-valued probability 
space {A,E). We assume: 

• All Xi 's have the same *- distribution with respect to E and their 
first moments vanish, i.e. E[Xi] = 0. 

• The Xi are *-free with amalgamation overM^iC) (which means 
that the *-algebras Xi, generated by Mm(C) and Xi, are free with 
respect to E). 

• We have sup < 00. 

Then the sequence {Sn)nm defined by 




n G N 
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converges to an operator-valued semicircular element s . Moreover, we 
can find k>0,c>1,C>0 and G N such that 

\GJb) - Gs (b)\\ < C^\\b\\ for all b e Q and n>N, 

'n 



where 

:= jfe G GL„(c) I iir^ii < \\b\ ■ liri < c} 

and where Gg and Gs„ denote the operator-valued Cauchy transforms 
of s and of Sn, respectively. 

Applying this operator- valued statement to our multivariate problem 
gives the following main result on a multivariate free Berry Esseen 
theorem. 

Theorem 1.2. Let (x-^^)^^j^, i G N, he free and identically distributed 
sets ofd self-adjoint non-zero random variables in some non- commutative 
C* -probability space {C,t), with t faithful, such that the conditions 

r(xf ^) =0 for k=l,...,d and alii eN 



and 



II (fc)ii ^ 
sup max \\xl || < oo 



k=l, 



are fulfilled. We denote by Ti = {<Jk,i)ki=i> where a^^i '■= T{x[''^xf^), 
their joint covariance matrix. Moreover, we put 

1 " 

Si^^ := ^ V ^ fork = l,...,d and all n G N. 

^ i=i 

Then {Si^\ . . . , Si""^) converges in distribution to a semicircular family 
(si, . . . , Sd) of covariance S. We can quantify the speed of convergence 
in the following way. Letp be a (not necessarily self-adjoint) polynomial 
in d non-commutating variables and put 

P„ := . . . , and P := p{s^, ...,s,). 

Then, there are constants C > 0, R > and N & N (depending on the 
polynomial) such that 

\Gp{z) - Gp {z)\ < G^ for all \z\ > R andn>N, 

\/n 



where Gp and Gp^ denote the scalar-valued Cauchy transform of P and 
of Pn, respectively. 

In order to deduce estimates for the Kolmogorov distance A(P'„, P) 
one has to transfer the estimate for the difference of the scalar-valued 
Cauchy transforms of P„ and P from near infinity to a neigborhood 
of the real axis. A partial solution to this problem was given in the 
appendix of [13], which we will recall in Section 4. But this leads to the 
still unsolved question, whether p{si, . . . ,Sd) has a continuous density. 
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We conjecture that the latter is true for any selfadjoint polynomial in 
free semicirculars, but at present we are not aware of a proof of that 
statement. 

The paper is organized as follows. In Section 2 we recall some basic 
facts about holomorphic functions on domains in Banach spaces. The 
tools to deal with matrix-valued Cauchy transform will be presented 
in Section 3. Section 4 is devoted to the proof of Theorem 11.11 and 
Theorem 11.21 

2. Holomorphic functions on domains in Banach spaces 

For reader's convenience, we briefly recall the definition of holomor- 
phic functions on domains in Banach spaces and we state the theorem 
of Earle-Hamilton, which will play a major role in the subsequent sec- 
tions. 

Definition 2.1. Let {X, \\ ■ \\x), {Y, \\ ■ be two complex Banach 
spaces and let D C X be an open subset of X. A function f : D ^ Y 
is called 

• strongly holomorphic, if for each x E D there exists a bounded 
linear mapping Df{x) : X Y such that 

^^ Jf{^ + y)-f{^)-Dfi^)y\\Y Q 

y^o \\y\\x 

• weakly holomorphic, if it is locally bounded and the mapping 

X^4>{f{x + Xy)) 

is holomorphic at A = for each x E y eY and all continu- 
ous linear functionals : F — ?■ C. 

An important theorem due to Dunford says, that a function on a do- 
main (i.e. an open and connected subset) in a Banach space is strongly 
holomorphic if and only if it is weakly holomorphic. Hence, we do not 
have to distinguish between both definitions. 

Definition 2.2. Let D he a, nonempty domain in a complex Banach 
space (X, II ■ II) and let / : D — )■ be a holomorphic function. We say, 
that f{D) lies strictly inside D, if there is some e > such that 

Be{f{x)) C D for all X G D 

holds, whereby we denote by Br{y) the open ball with radius r around 

y- 

The remarkable fact, that strict holomorphic mappings are strict 
contractions in the so-called Caratheodory-Riffen-Finsler metric, leads 
to the following theorem of Earle-Hamilton (cf. [3j), which can be seen 
as a holomorphic version of Banach's contraction mapping theorem. 
For a proof of this theorem and variations of the statement we refer to 
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Theorem 2.3 (Earle-Hamilton, 1970). Let ij) D C X be a domain in 
a Banach space {X, \\ ■ ||) and let f : D ^ D be a bounded holomorphic 
function. If f{D) lies strictly inside D, then f has a unique fixed point 
in D. 

3. Matrix-valued spectra and Cauchy transforms 

The statement of the foUowing lemma is weU-known and quite sim- 
ple. But since it turns out to be extremely helpful, it is convenient to 
recall it here. 

Lemma 3.1. Let || • ||) be a complex Banach- algebra with unit 1. 
If X & A is invertible and y & A satisfies \\x — y\\ < c-p^ for some 
< a < 1, then y is invertible as well and we have 

M IM 1 n IM 

1 — a 

Proof We can easily check that 

oo 

J2{x-\x-y)rx-^ 

n=0 

is absolutely convergent in A and gives the inverse element of y. More- 
over we get 

oo ^ 
n=0 

which proves the stated estimate. □ 

Let (C,r) be a non-commutative C*-probability space, i.e., C is a 
unital C*-algebra and r is a unital state (positive linear functional) 
on C; we will always assume that r is faithful. For fixed m G N we 
define the operator- valued C*-probabihty space A := M„(C) with 
conditional expectation 

E -.^idm^r : ^^M^(C), 6 c i-^ r(c)6, 

where we denote by Mm(C) the C*-algebra of all m x m matrices 
over the complex numbers C. Under the canonical identification of 
Mm(C) ® C with Mjn{C) (matrices with entries in C), the expectation 
E corresponds to applying the state r entrywise in a matrix. We will 
also identify b G M„j(C) with 6 (g) 1 G ^. 

Definition 3.2. For a E A = M^(C) we define the matrix-valued 
resolvent set 

Pm{,o) := {h G M„j(C) I 6 — a is invertible in A} 

and the matrix-valued spectrum 

am{a) := M„(C)\p„(a). 
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Since the set GL(^) of all invertible elements in A is an open subset 
of A (cf. Lemma [3 -ip . the continuity of the mapping 

fa : Mm{C)^A, b^b-a 

implies, that the matrix-valued resolvent set Pm{o) = fa^i^^i.-^)) 
element a G ^ is an open subset of Mm(C). Hence, the matrix- valued 
spectrum (Jm{o) is always closed. 

Although the behavior of this matrix-valued generalizations of the 
classical resolvent set and spectrum seems to be quite similar to the 
classical case (which is of course included in our definition for m = 1), 
the matrix valued spectrum is in general not bounded and hence not a 
compact subset of Mm(C). For example, we have for all A G C, that 

a^(Al) = {be M„(C) I A G aM^colf*)}, 

i.e. crm(Al) consists of all matrices b G Mm(C) for which A belongs to 
the spectrum crM„(c)(^)- Particularly, a"m(Al) is unbounded for m > 2. 

In the following, we denote by GLm(C) := GL(Mm(C)) the set of all 
invertible matrices in Mm(C). 

Lemma 3.3. Let a E A be given. Then for all b G GLm(C) the 
following inclusion holds: 

{A6 I A G PA{b'^a)} ^ Pm{a) 

Proof. Let A G p_4(6^^a) be given. By definition of the usual resolvent 
set this means that Al — b~^a is invertible in A. It follows, that 

\b- a = b{Xl-b-^a) 

is invertible as well, and we get, as desired, Xb G Pm{o)- D 

Lemma 3.4. For all ^ a E A we have 

{6GGL,„(C)|||ri<^}cp^(a) 

and 

aM n GL„(C) C jfe G GL^(C) | \\b-'\\ > 

Proof. Obviously, the second inclusion is a direct consequence of the 
first. Hence, it suffices to show the first statement. 

Let b G GLm,(C) with \\b^^\\ < ^ be given. It follows, that h := 
1 — b^^a is invertible, because 

||1 - h\\ = \\b-^a\\ < \\b'^\\ ■ \\a\\ < 1. 

Therefore, we can deduce, that also 

(3.1) fe-a = 6(l-r^a) 

is invertible, i.e. b G Pm{o)- This proves the assertion. □ 
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The main reason to consider matrix-valued resolvent sets is, that 
they are the natural domains for matrix-valued Cauchy transforms, 
which we will define now. 

Definition 3.5. For a G ^ we call 

Ga-. p™(a)^M„(C), b^E[{b-a)-'] 
the matrix-valued Cauchy transform of a. 

Note that Ga is a continuous function (and hence locally bounded) 
and induces for all bo G Pm{o)i b G Mm(C) and bounded linear func- 
tional ^ : v4 — )■ C a function 

which is holomorphic in a neighborhood of A = 0. Hence, Ga is weakly 
holomorphic and therefore (as we have seen in the previous section) 
strongly holomorphic as well. 

Because the structure of Pm{o) and therefore the behavior of Ga 
might in general be quite complicated, we restrict our attention to a 
suitable restriction of Ga- In this way, we will get some additional 
properties of Ga- 

The first restriction enables us to control the norm of the matrix- 
valued Cauchy transform on a sufficiently nice subset of the matrix- 
valued resolvent set. 

Lemma 3.6. Let ^ a & A be given. For < 9 < 1 the matrix valued 
Cauchy transform Ga induces a mapping 

9 1 



9 Hal 



Ga-. {be GL„(C) I liri < ^-^l ^ jft G M^(C) I < 
Proof. Lemma 13.41 (c) tells us, that the open set 

u--= {be GL„(c) I ||ri < 9 ■ ^1 

is contained in Pm{a), i.e. Ga is well-defined on U . Moreover, we get 
from 

00 

ib - a)-i = (1 - b-^ayh-^ = [b-^afb-^ 

n=0 

and hence 

00 ^ ^ 

(3.2) ||G.(6)|| < \\{b-a)-'\\ < lir^ll J](||f>-i||||a||)" < 



1-9 \\a\ 

n=0 " ' 



for all b eU . This proves the claim. □ 

To ensure, that the range of Ga is contained in GLm,(C), we have to 
shrink the domain again. 
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Lemma 3.7. Let ^ a E A he given. For < 9 < 1 and c > 1 we 

define 

:= {fee GL,„(C) I \\h-^\\ < ^ - ll^ll • ll^^i < 

and 

Q' :=|bGGL„(C)|||b||< 
// the condition 



1-e 



a 

< - 



1-9 c 

is satisfied for some < a < 1, then the matrix-valued Cauchy trans- 
form Ga induces a mapping : i7 — )■ and we have the estimates 

(3.3) < \\{b-a)-'\\<-^--^ for all ben 

1 — 9 \\a\\ 

and 

(3.4) \\Ga(b)-^\\ < ^— -m for all ben. 

1 — a 

Proof. For all 6 G ^2 we have 

oo 

G,{b) - b-^ = E[{b - a)-' - b-'] =E^J2 [b^'cLyb-' 
which enables us to deduce 



n=l 



l|G.(6)-r'|| < ||6-'|| Y, < Y^gWo-'W < < - 

Using Lemma 13.1^ this implies Ga{b) e GLm(C) and fl3.4p . Since we 
already know from fl3.2p in Lemma I3.6[ that f l3.3p holds, it follows 
Ga{b) e n' and the proof is complete. □ 

Remark 3.8. Since domains of our holomorphic functions should be 
connected it is necessary to note, that for k > and c > 1 

n = {be GLmic) I ||ri < K, \\b\\ ■ \\b-^\\ < c] 

and for r > 

n' = {beGLm{C) I \\b\\<r} 

are pathwise connected subsets of Mm{C). Indeed, if 61,62 £ GLm(C) 
are given, we consider their polar decomposition bi = UiPi and 62 = 
U2P2 with unitary matrices Ui,U2 e GLm(C) and positive-definite Her- 
mitian matrices Pi, P2 e GLm(C) and define (using functional calculus 
for normal elements in the C*-algebra Mm(C)) 

7 : [0, 1] ^ GL„(C), t ^ Ul-'Pl-'UlPl 

Then 7 fulfills 7(0) = 61 and 7(1) = 62, and 7([0, 1]) is contained in n 
and n' if 61, 62 are elements of n and n' , respectively. 
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Since the matrix-valued Cauchy transform is a solution of a special 
equation (cf. [T2]), we will be interested in the following situation: 

Corollary 3.9. Let rj : GLm(C) — > Mm(C) be a holomorphic function 
satisfying 

\\t]{w)\\ < M\\w\\ for all w G GL^(C) 
for some M > 0. Moreover, we assume that 

hGa{h) = l + r]{Ga{h))Ga{h) forallben 
holds. Let < 9, a < 1 and c > 1 be given with 



rl Hal 

< crmm 



{c' "m } 



1 - 

and let Q and Q' be as in Lemma \3. 71 
Then, for fixed b ^ Q, the equation 

(3.5) bw = 1 + ri{w)w, w & Q' 

has a unique solution, which is given by w = Ga{b). 

Proof. Let b E Q he given. For all w E Q' we get 

||r7H||<M||«;||<^-^ 

and therefore 

lir^Hii < liriiir^Hii < ■ ft^ ■ ^ < < i- 

1 — U ||a||^ 

This means, that 1 — b^^rjiyo) and hence b — ri{w) is invertible with 
\\{b-ri{w))-'\\<\\b-^\\{l-b-'ri{w)r'\\ 

oo 
n.=0 

e 1 



< 



l-Oa \\a\y 

and shows, that we have a well-defined and holomorphic mapping 
J- : fi' ^ M^(C), w^ib- r]{w))-^ 

with 

II^HII = 11(6- r7H)-i < • A < • A 

\ — da \\a\\ 1 — u \\a\\ 

and therefore J^{w) G Vt' . 

Now, we want to show that J-'{Vt') lies strictly inside Vt' . We put 

1 - ^ \ 6 1 
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and consider w & Q' and u G Mm(C) with \\u — J^{w)\\ < e. At first, 
we get 

\\b - v{w)\\ < \\b\\ + Mw)\\ < ll&ll + o ■ ^ ll^ll + 

||a|| i — (7 

and thus 

||m-(&-?7(-u;))~"^|| = \\u-J^{w)\\ < e < ^- ^ <^ ^ 



2 \\b\\ + a\\a\\ " 2 \\b - 7]{w) 
which shows u G GLm(C), and secondly 

ll^ll = \\u-{b-r]{w))-'\\ + \\J^{w)\\ 

1-e e 1 



< 



l-Oa i-e 
e 1 



9 \\a\ 



which shows u eVL' . 

Let now w he a, solution of (13. 5p . This implies that 

w^^J^{w) = w^^{b — r]{w))~^ = {bw — 'r]{w)w) ""^ = 1, 

and hence T{w) = w. Since : Q' Q' is holomorphic on the 
domain Q' and J^{Q') lies strictly inside Q', it follows by the Theorem 
of Earle-Hamilton, Theorem 12. 3[ that J-" has exactly one fixed point. 
Because Ga{b) (which is an element of Q' by Lemma [3]Tj) solves (13. 5p by 
assumption and hence is already a fixed point of J-", it follows w = Ga{b) 
and we are done. □ 

Remark 3.10. Let [A! , E') be an arbitrary operator- valued C*-probability 
space with conditional expectation E' : A! ^ Mm(C). This provides 
us with a unital (and continuous) ^-embedding l : Mm(C) — )■ A! . In 
this section, we only considered the special embedding 

L : M^(C) ^ A, b^b®l, 

which is given by the special structure A = Mm{C) ® C. But we can 
define matrix- valued resolvent sets, spectra and Cauchy transforms also 
in this more general framework. To be more precise, we put for all 
aeA' 

Pm{(i) '■= {b G Mm(C) I L{b) — a is invertible in A'} 

and o"m(a) := Mm(C)\pm(a) and 

Ga-. p„(a) ->M„(C), b^E'[{L{b)-a)-']. 

We note, that all the results of this section stay valid in this general 
situation. 
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4. Multivariate free central limit theorem 

4.1. Setting and first observations. Let (Xj)jgN be a sequence in 
the operator-valued probability space {A, E) with A = M^iC) = 
Mm(C) ® C and E = id®r, as defined in the previous section. We 
assume: 

• All XiS have the same *-distribution with respect to E and 
their first moments vanish, i.e. = 0. 

• The Xi are *-free with amalgamation over Mm(C) (which means 
that the *-algebras Xi, generated by Mm(C) and Xi, are free 
with respect to E). 

• We have sup \\Xi\\ < oo. 

If we define the linear (and hence holomorphic) mapping 

T] : M^(C) -> M^(C), b ^ E[X,bX^], 
we easily get from the continuity of E, that 

\\r]{b)\\ < f sup Vll for all b G M^(C) 

holds. Hence we can find M > such that ||?7(^))|| < holds for all 

b G Mm(C). Moreover, we have for all G N and all 6i, . . . , 6^ G Mm(C) 

/ N fc + l 

sup \\E[XibiXi . . .bkXi]\\ < I sup ) ■ • • \\bk\\- 

Since (Xj)jgN is a sequence of centered free non-commutative random 
variables. Theorem 8.4 in [15j tells us that the sequence {Sn)nm defined 
by 



1 

•= :7^IZ^'' 



^^ G N 



n 

i=l 



converges to an operator- valued semicircular element s. Moreover, 
we know from Theorem 4.2.4 in [T2] that the operator-valued Cauchy 
transform Gs satisfies 

bGs{b) = 1 + r]{Gs{b))Gs{b) for all b G Ur, 

where we put Ur := {b G GLm(C) | \\b^^\\ < r} C pm{s) for all suitably 
small r > 0. 

By Proposition 7.1 in [S], the boundedness of the sequence (Xj)jgp^ 
guarantees boundedness of {Sn)n&i as well. In order to get estimates 
for the difference between the Cauchy transforms Gg and Gs„ we will 
also need the fact, that {Sn)nm is bounded away from 0. The precise 
statement is part of the following lemma, which also includes a similar 
statement for 

1 1 " 

:= Sn ^Xi = ^ V X,- for all n G N and 1 < « < n. 
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Lemma 4.1. In the situation described above, we have for all n E N 
and all 1 < i < n 



\\Sn\\>\\a\\"^ and ||5W||>Wl-- 

V n 

where a := E[X*Xi\ e M^(C). 

Proof. By the *-freeness of Xi, X2, . . . , we have 

E[X*X,] = E[X*] . E[X,] = 0, for t + j 



and thus 

WSJ' = \\s:sn\\ > \\E[s:sr, 

Similarly 



>\\E[(s^^rsii^]\\ 



E[SnSn\ E[X*Xi 

n 



n- 1 



\a\ 



n 



□ 



-1 



which proves the statement. 
We define for n G N 

Rn •■ Pm{Sn) ^ A, b ^ {b - 5„)"^ 

and for 77. e N and 1 < i < n 

Lemma 4.2. For all n eN and 1 < i < n we have 

(4.1) R^{b) = i?|:l(6) + + -R^{b)X,R^^{b)X,R^^{b) 

and 
(4.2) 



R^{b) = i?[:l(6) + -^R^^{b)X,Rr,ib) 



/or a//6ep^(5„)np^(5[:l). 
Proof. We have 
(6-5„)/246)(6-S'W) = 



= (6-5„) + -^(6-5W)i2W(6)X, 

= (6 - 5„) + ^(6 - S^)R^^{b)X, + -X,R§{b)X, 

\/ Th Til 
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which leads, by multiphcation with Rn{b) = {b — Sn)~^ from the left 
and with Rn\b) = {b — Sn^)~'^ from the right, to equation (14. ip . 
Moreover, we have 

{b - (6 - 5„) = 6 - = {b- Sr.) + 



which leads, by multiplication with Rn{b) = {b — Sn) ^ from the right 
and with = {b - SH^)'^ from the left, to equation (H^D- □ 

Obviously, we have 

:= Gs„ = Eo i?„ and := G^,, =EoRI^. 

4.2. Proof of the main theorem. During this subsection, let < 
^, o" < 1 and c > 1 be given, such that 

(4.3) < crmm <^ -, 



i-e Ic' M 

holds. For all n G N we define 

« • ■ 1 1 1 
Kn '■= c/mm 

and 



sr\\Sn\y\\sW"'^\\s^\ 



fi„ := {b G GL^(C) I liri < • \\b-^\\ < c}. 

Lemma [3.41 shows, that fi„ is a subset of Pm{Sn)- 

Theorem 4.3. For all 2 < n & N the function Gn satisfies the follow- 
ing equation 

An{b)Gnib) = 1 + ?7(G„(6))G„(fe), b e n^, 

where 

An : Qn^MmiC), b^b- Qn{b)Gn{b)-\ 
with a holomorphic function 

e„: 1^„^M„(C) 

satisfying 

sup ||e„(6)|| < 



b&n„ y'n 
with a constant G > 0, independent of n. 

Proof, (i) Let n G N and b G Pm{Sn) be given. Then we have 
SMb) = bRn{b) -{b- Sn)Rn{b) = bRnib) - 1 

and hence 

E[SnRn{b)] = E[bRn{b) - l] = 6G„(6) - I. 
(ii) Let n G N be given. For all 



b G pm,„ := pm{Sn) n p|p™(5'i*l 



i=l 
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we deduce from the formula in (14. ip . that 

1 " 

+ -E[XiRr,{b)XiR^^{b)X,R^^{b)] 
= -f2(E[X.Rl^m..R^^ib)] + -^E[X,R^ib)X^R^m^R^ 

i=l \ V / 

1 " 



n 

i=l 



where 

: pUSn)npUSli^) ^M^C), b^ ^E[X,R^{b)X,R^^{b)X,R^^{b)], 



n 



There we used the fact, that, since the (Xj)jgp^ are free with respect to 
E, also Xi is free from Rn\ and thus we have 

E[X,R^:}{b)]=E[X,]E[R^^{b)]=0 

and 

E[X,R^:^{b)X,R^:^{b)] = E[X,E[R^^{b)]X,]E[Rf}{b)]. 
(iii) Taking (14. 2 p into account, we get for all n G N and 1 < i < n 

Gn{b) = E[R^{b)\ = E[R}^:^{b)\+^E[R^^{b)X,Rr.{b)\ = G^:}{b)-r^^^,{b) 



n 



and therefore 

G^:^{b) = Gn{b)+r%{b) 
for all b G Pm{Sn) H PmiSn), where we put 

rZ2- PmiSn)n pUSH^)^ MUC), b^-^E[R^^ib)X,R^ 
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(iv) The formula in (iii) enables us to replace Cl^' in (ii) by Gn- Indeed, 
we get 



E[SMb)] = + 

i=l 
i=l 

1 



n 

i=l 



for all b G pm,ni where the function 
is defined by 

r%{b) ■■= viG^mrlM + il{r%{h))G^{h) + v{r%{h))r%{b) + r^M- 
(v) Combining the results from (i) and (iv), it follows 

bGn(b) - 1 = E[SnRn{b)] = r^{Gn(b))Gn{b) + Qnib), 
where we define 

1 " 
n ^ — ^ 

i=l 

Due to (14. 3p , Lemma 13.41 and Lemma 13.71 show that C „ and 
Gn{b) e GL„(C) for b e Qn- This gives 

(f>-e„(b)G'„(6)-i)G'„(6) = l + v{Gn{b))Gr.{b) 

and hence, as desired, for all 6 G fi„ 

An(6)G'„(6) = l + r7(G„(6))G„(6). 

(v) The definition of Qn gives, by Lemma 13.31 and by Lemma 14. 1[ the 
following estimates 

\\G4b)\\ < \\Rnm < ■ ^ < ■ b G a„ 

and 

Therefore, we have for all 5 G fi„ by (ii) 

iirWi(6)ii < -^wx^rimbMR^Mr < ^^(t^^) V.ir 
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and by (iii) 

and finally by (iv) 

11413(6)11 < 2M||G„(6) II 1141,(6) II + M||ri'y6)|p + 11^,(6)11 

< — 

In 



for all h G where C > is a constant, which is independent of n. 
Hence, it follows from (v) that 

sup ||e„(6)||<-^. 

□ 

The definition of n„ ensures, that 

G := a : ^ M^(C) 

satisfies 

hG{h) = 1 + r](G(b))G(b) for all b e Q, 

where 

Q := {& e GL^(C) I \\b-'\\ < ^ ■ ll^ll ■ ll^^i < c} - 
We choose 

(4.4) < 7 < and < r < (1 - 7)^ 

(note, that < 7 < 1) and we put c* := c — (1 + 0)7, which fulfills 
clearly 1 < c* < c. Since we have 9* < 9 and c* < c, we see 

9* , 9 

c < < a 



1-9* 1-9 
and hence 

Finally, we define 

1 1 



K* ■=9* min 



UNI'||5n||'||5W||'---'||^W||'' 
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and 

:= [b e GL„(C) I iirl < <, \\b\\ ■ liri < c*} c n^,. 

Corollary 4.4. There exists G N such that 

AniK) ^ for all n>N. 



Proof. Since we have by Theorem 14.31 

C 



sup ||0„(6)|| < 

for all 2 < n G N, we can choose an G N such that 

sup ||e„(6)||<^(l-a) 

holds for all n> N. Now, we get for all & G fi* : 

(i) An{b) is invertible: Since (13. 4p gives 

\\GrXb)-'\\ < -r^m forall6Gfi„, 
i — O" 

we immediately get 
\\An{b) - fell < ||e.(6)||||G„(b)-i < 7^ < < 

(ii) We have ||A„(f>)^^|| < /c„: Using Lemma l3.1[ we get from (i) 
that 

l|A„(fe)-^||< ^116-^11 <T^<«:„. 
1 — 7 1 — 7 

(iii) We have ||A„(fe)||||A„(6)-i|| < c: Using 

l|AnW-fe|| <7^ 

from (i) and 

\\K{b)-'\\ < -r^\\b-'\\ 
1 - 7 

from (ii), we get 
||A„(6)||||A„(r^|| < (||fe|| + ||A„(6)-6||)||A„(6)-^|| 
< (, , 7^ 1 



c* / 1 — 7 

c* +7 
< < c. 

1-7 

Finally, this shows A„(6) G i7„,. □ 
Corollary 4.5. For all n > N we have 

Gn{b) = G{An{b)) for all b e Q*^. 
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Proof. For all ?7, G N we define 

n'^:= |6gGL„(C)|||6||<^}. 

Let n > N and G ^7* be given. We know, that 

A„(6)G'(A„(b)) = 1 +r7(G(A„(fo)))G(A„(6)) 
holds, i.e. w = G'(A„(6)) G Q'^ ^ solution of the equation 

An{b)w = 1 + 'r]{w)w, w G Vl'^. 
Combining (14. 3 p with Lemma [4. ![ we get 

< a mm < -, — — - ) < a mm < -, — — — : n G N 

i-e \c M ) - \c M 

Hence, the equation above has, by Theorem 13. 9[ the unique solution 
w = Gn{b) G Q'n- This implies, as desired, Gn{b) = G'(A„(6)). □ 

Corollary 4.6. For all n > N we have 

\\Gib) - G„(6)|| < C'^\\b\\ for all b G fi^, 



where C > is a constant independent of n. 

Proof. For all 6 G fi* C n„ C we have 
G{b)-Gn{b) = Gib) - G{kn{b)) 

= E[{b - s)-' - {K(b) - s)-'] 
= E[{b-.s)~\K{b)-b){K{b) 

and therefore by (13.41) . which gives 

\\Gn{b)-'\\ < -^\\b\\ for all 6 Gl^:, 
i — a 

and (since A„(6) G fi„ C by 

\\G{b) - G^m < ||(b-s)-i||.||A„(6)-6||-||(A„(fo)-s) 

9 1 \2 



- (t^'m) ■ \\^r.m ■ \\Gn{b) 



-1 1 



'n ' 



where 

This proves the corollary. □ 

We recall, that the sequence (Xj)jgN is bounded, which implies bound- 
edness of the sequence {Sn)neN as well. This has the important conse- 
quence, that 

111 1 
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for some k* > 0. If we define 

Q* ■= jfe G GL„(c) I ||ri < K*, \\b\\ ■ \\b-^\\ < c*y 

we easily see fl* C f2* for all n G N. Hence, by renaming Q* to Q etc., 
we have shown our main Theorem 11.11 

We conclude this section with the following remark about the geo- 
metric structure of subsets of Mm(C) like Q. 

Lemma 4.7. For k > and c > 1 we consider 

:= |6 G GL„(C) I lirl < K, \\b\\ ■ \\b-^\\ < c}. 

For A, /i G C\{0} we define 



G GL, 



If j: < holds, we have A(A, fi) & Q for all 

(4.6) max | — , —1 < lAI < c|/i|. 

Ik c J 

Particularly, we have for all |A| > j:, that Al G fi. 

Proof. Let fi G C\{0} with - < be given. For all A G C\{0}, which 
satisfy ( 14. 6p . we get 

||A(A,/x)-i = ||A(A-\/i-i)|| = max{|Ar\ \fi\-'} < k. 

and 

||A(A,/i)|| • ||A(A,;u)^"'"|| = max I |A|, } ■ max I |A|^\ 





1^ 


. 


• °\ 







/i . 


. 


A(A,/i) := 








Vo 


. 





fl/^l|Ar\ 


if 


|A| 


< 


1^1 


liAii^r\ 


if 


|A| 


> 





< c, 

which implies A(A,//) G fi. In particular, for A G C\{0} with |A| > - 
we see that /i = A fulfills (gj]) and it follows Al = A(A, A) G fi. □ 

4.3. Application to multivariate situation. 

4.3.1. Multivariate free central limit theorem. Let {xf^)'^^'^^ z G N, be 
free and identically distributed sets of d self-adjoint non-zero random 
variables in some non-commutative C*-probability space (C, r), with r 
faithful, such that 

r(xf^) = for fc = l,...,dand alH G N 
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and 



(4.7) sup max \\x\''^\\ < oo. 

jgN k=l,...,d 

We denote by S = {<Tk,i)ii=i, where ak^i := T{xf'^xf^), their joint co- 
variance matrix. Moreover, we put 

1 " 

Si^^ -=^^2 ^i''^ for A; = 1, . . . , d and all ri G N. 



'n 

1=1 

We know (cf. [II]), that {Si^\ ...,SL^^) converges in distribution as 
?7, — > oo to a semicircular family (si, . . . , Sd) of covariance S. For no- 
tational convenience we will assume that Si, . . . , live also in (C, r); 
this can always be achieved by enlarging (C,r). 

Using Proposition 2.1 and Proposition 2.3 in for each polynomial 
p of degree g in d non-commuting variables, we can find m G N and 
ai, . . . , ttrf G Mm(C) such that 

Xl-piS^^\...,Si'^) and Al-p(si,...,s,) 

are invertible in C if and only if 

A(A,1)-S'„ and A(A, 1) - s, 

respectively, are invertible in ^ = Mm(C). The matrices A(A, 1) G 
Mm{C) were defined in Lemma HTTl and Sn and s are defined as follows: 

d 

:= ^ flfc ® Si^^ G A for all n G N 



fc=i 



and 



If we also put 



s := ^ Uk 1^ Sk e A. 



k=l 



Xi := ^ Cfc (g) xf^ G ^ for all i G N, 

fc=i 

then we have 

1 " 

'S'n = —;= > ^i. 



n 

■ 1 



We note, that the sequence (Xj)jgpj is *-free with respect to the con- 
ditional expectation E : A = Mm(C) — )■ Mm(C) and that all the Xj's 
have the same *-distribution with respect to E and that they satisfy 
= 0. In addition, (14. 7p implies supj^^j ||Xj|| < oo. Hence, the 
conditions of Theorem II. II are fulfilled. But before we apply it, we note 
that {Sn)n&i converges in distribution (with respect to E) to s, which 
is an Mm(C)-valued semicircular element with covariance mapping 

f]-. M^(C) ^M„(C), h^E[shsl 
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which is given by 

d d 

r]{b) = E[sbs] = ^ id ®r[(afc (g) l)(az ® si)] = ^ akhaiak,i- 

k,l=l k,l=l 

Now, we get from Theorem 11.11 constants > 0, c* > and C" > 
and iV G N such that we have for the difference of the operator-valued 
Cauchy transforms 

Gs{b) ■.= E[{b-s)-'] and Gs^b) := E[{b - Sn)-'] 
the estimate 

\Gs(b) - GsJb)\\ < C"^||6|| for all beQ* andn> N, 

'n 



where we put 

:= \h e GL„(C) I \\b-^\\ < k\ \\b\\ ■ \\b-^\\ < c*y 

Moreover, Proposition 2.3 in j5] tells us 

(Al - p{Sl^\ Si'-)))-' = (tt ® idc) ((A(A, 1) - S^)-') 

and 

(Al - p{si, Sd)) = ® idcO ((A(A, 1) - s)-') , 

where tt : Mm(C) — )■ C is the mapping given by 7r((ajj)jj=i^...^m) := di,!- 
Since to [n^idc) = tto E, this implies a direct connection between the 
operator-valued Cauchy transforms of Sn and s and the scalar-valued 
Cauchy transforms of P„ := p{Sn \ . . . , Sn^) and P := p{si, . . . ,Sd), 
respectively. To be more precise, we get 

GpJA) := r[(A - P„)-i] = 7r(G5„(A(A, 1))) 

and 

Gp{X) ■.= T[{X-P)-'] = n{Gs{A{X,l))) 

for all A G Pc{Pn) and A G Pc{P), respectively. 

If we choose /i G C such that |/i| > 4r holds, it follows from Lemma 
K7\ that A{X,fi) G n* is fulfilled for all A G where A{fi) C C 

denotes the open set of all A G C satisfying ( 14. 6p . i.e. 

A{fi) := |a G C I max ^} < |A| < c*|/i|}. 

If we apply Proposition 2.1 and Proposition 2.2 in |5J to the polyno- 
mial -p^p (which corresponds to the operators j^Sn and j^S), we easily 
deduce that 

Al - -^^p{S^^\ Si'^) and Al - -^,p{su 

are invertible in C if and only if 

A{X, n) — Sn and A{X, p) — Sn, 
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respectively, are invertible in A. Moreover, we have 
/i^-^Gp„(A/i^"i) = 7r(G5„(A(A,/x))) 

and 

p'~'GpiXfi'-') = n{GsiAiX,p))) 
for all A e pcijJ^Pn) and A G pdjJ^P), respectively. 

Particularly, for all A G A{p) we get A(A, p) G Q* and hence A G 
PcijJ^Pn) n PcijJ^P) for all n> N. Therefore, Theorem 11.11 implies 

\pr'\Gp{Xp^-') - Gp„{Xp^-')\ = \n{Gs{A{X,p))-GsAHKp)))\ 

< \\G,{A{X,p))-Gs,AH\p))\\ 

< C'^\\A{X,p)\\ 



< C"^max{|A|,|^|} 

'n 



< C'c*\X' ^ 



n 



\Gp{Xp^-') - Gp,XXp<^~')\ < G'c*^\Xp^-'\ 



and hence 



This means, that 

Gp{z)-Gp„{z)\<G'c*^\z\ 

n 



holds for all 2 G C with G A{p) and all n > N. By definition of 
A{p), we particularly get 

\Gp{z) - Gp„{z)\ < G^ for all — < \z\ < c*\p\^ and n > A^, 

'n c* 



where we put G := G'{c*)'^\p\3 > 0. Since z i-> Gp{z) - Gp^{z) is 
holomorphic on {2; G C | \z\ > R} for R := > and extends 

holomorphically to 00, the maximum modulus principle gives 

\Gp{z) - Gp^{z)\ < G^ for all \z\ > Randn> N. 



n 



This shows Theorem [L21 



4.3.2. Estimates in terms of the Kolmogorov distance. In the classical 
case, estimates between scalar-valued Cauchy transfoms can be estab- 
lished (for self-adjoint operators) in all of the upper complex plane and 
lead then to estimates in terms of the Kolmogorov distance. In the 
case treated above, we have a statement about the behavior of the dif- 
ference between two Cauchy transforms only near infinity. Even in the 
case, where our operators are self-adjoint, we still have to transport 
estimates from infinity to the real line, and hence we can not apply 
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the results of Bai |Jj directly. A partial solution to this problem was 
given in the appendix of [13] with the following theorem, formulated in 
terms of probability measures instead of operators. There we use the 
notation for the Cauchy transform of the measure /i, and put 

D^:={zeC \ lm{z) > 0, \z\ > R}. 

Theorem 4.8. Let fi be a probability measure with compact support 
contained in an interval [—A, A] such that the cumulative distribution 
function J-"^ satisfies 

\J^,i{x + t) - Tf,{x)\ < p\t\ forallx,teR 

for some constant p > 0. Then for all R > and j3 G (0, 1) we can 
find > and nio > such that for any probability measure u with 
compact support contained in [—A, A], which satisfies 

sup \G^{z)-G,{z)\<e-"' 

for some m > rriQ, the Kolmogorov distance A{fi,u) := snp \J^n{x) — J^^{x)\ 
fulfills 

A(/i,//)<e^. 

mP 

Obviously, this leads to the following questions: First, the stated 
estimate for the speed of convergence in terms of the Kolmogorov dis- 
tance is far from the expected one. We hope to improve this result in 
a future work. Furthermore, in order to apply this theorem, we have 
to ensure that p{si, . . . ,80) has a continuous density. As mentioned in 
the introduction, it is a still unsolved problem, whether this is always 
true for any self-adjoint polynomials p. 
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